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Chapter  1 
Introduction 


This  report  will  describe  two  improvements  to  the  “Electromagnetic  Surface 
Patch"  or  ESP  code.  ESP  is  a  user  oriented  computer  code,  baaed  upon  the 
method  of  moments  (MM)  (1),  for  the  analysis  of  thin  wires  and  polygonal 
plates.  ESP  has  been  distributed  to  over  100  users  in  industry,  government, 
and  universities. 

The  latest  version,  ESP  Version  III,  was  released  in  May  1987  (2j.  In 
brief,  ESP  Version  III  can  treat  geometries  consisting  of: 

1.  thiu  wires  with  finite  or  perfect  conductivity  and  lumped  loads, 

2.  perfectly  conducting  polygonal  plates, 

3.  wire/ plate  junctions  (at  least  0.1  A  from  a  plate  edge) 

4.  plate/plate  junctions,  iucludiug  several  plates  which  intersect  along  a 
common  edge, 

5.  excitation  by  either  a  plane  wave  or  a  delta  gap  voltage  generator. 

ESP  III  is  capable  of  computing  most  of  the  quantities  of  engineering  in¬ 
terest  including: 

1.  current  distributions, 

2.  input  impedance,  input  admittance,  and  radiation  efficiency, 

3.  mutual  coupling, 
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4.  fax  zone  gain  patterns  (both  polarizations) 

5.  far  zone  plane  wave  scattering  patterns  (full  scattering  matrix). 

Notice  that  ESP  III  can  only  treat  perfectly  conducting  plates.  However, 
using  the  sheet  impedance  approximation  we  have  modified  the  code  so  that 
it  can  treat  electrically  thin  dielectric  plates.  The  dielectric  plates  may  be 
lossy.  Also  ESP  III  can  only  compute  far  zone  radiation  and  scattering 
patterns.  However  ESP  has  now  been  modified  to  be  able  to  compute  near 
zone  radiation  and  scattering  patterns,  including  the  radial  component  of 
the  field.  By  near  zone  patterns  we  mean  that  the  field  point  can  be  the  in 
near  zone  of  the  radiating  antenna  or  scatterer.  These  two  modifications 
have  been  incorporated  into  the  ESP  code,  and  will  become  part  of  ESP 
Version  IV.  All  of  the  numerical  data  in  this  report  was  run  with  this 
preliminary  version  of  ESP  IV.  Although  it  will  not  be  described  in  this 
report,  ESP  has  also  been  modified  so  that  it  can  efficiently  generate  data 
over  a  wide  frequency  bandwidth  {3j. 

In  this  report  all  fields  and  currents  are  time  harmonic,  with  the 
time  dependence  suppressed. 


2 


Chapter  2 

Modifications  to  ESP 


2.1  Near  Zone  Fields 

As  mentioned  above,  previous  versions  of  the  ESP  code  are  limited  to  com¬ 
puting  far  zone  patterns.  That  is,  the  field  point  is  at  an  infinite  radius 
from  the  radiating  or  scattering  body.  However,  we  have  now  increased  the 
generality  of  the  ESP  code  by  allowing  the  field  point  to  be  at  a  finite  ra¬ 
dius  from  the  body.  While  in  the  far  zone  only  the  transverse  components 
of  the  fields  are  nonzero,  in  the  near  zone  the  nonzero  radial  component  is 
computed. 

Cain  and  radar  cross-section  (RCS)  at  a  near  zone  field  point  will  now  be 
defined  so  that  they  arc  consistent  with  the  usual  far  zone  definitions.  We 
will  start  with  the  definition  of  gain.  Thus  we  consider  an  antenna  excited 
by  a  generator  and  which  radiates  an  electric  field,  E(r,0,  $),  which  is  a 
function  of  the  spherical  (r,  B,  d>)  coordinates  of  the  field  point.  Let  f?(r,  B ,  $} 
denote  a  given  polarization  of  the  vector  electric  field.  It  is  well  known  that 
in  the  far  zone  (i.e.,  in  the  limit  as  r  -•  oo)  the  radial  dependence  of  the 
field  is  simply 

«“*7  r. 

For  this  reason,  in  the  far  sour,  the  electric  field  is  usually  written  as 

£M,0)  =  —£>(«,*),  (2.1) 

f 

where  Ep  is  referred  to  as  the  far  zone  electric  field  and  is  independent  of 
r.  Front  Equation  2.1,  the  far  zone  electric  field  is  related  to  the  actual 
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electric  field  by 


EF(9t4>)  =  hmre>'£(r,M).  (2.2) 

The  far  zone  zone  directivity  of  an  antenna  in  a  given  direction  is  the 
ratio  of  the  radiation  intensity  in  that  direction  to  the  average  radiation 
intensity  of  the  antenna.  The  far  zone  gain  is  the  product  of  the  directivity 
and  the  radiation  efficiency.  The  radiation  efficiency  is  the  ratio  of  the  total 
power  radiated  by  the  antenna  to  the  total  power  input  to  the  antenna.  The 
far  zone  gain  is  simply  related  to  the  far  zone  electric  field  by 


Gp(9,d) 


!£WM)L* 

30Pm  ’ 


(2.3) 


where  Pxn  is  the  total  power  input  to  the  antenna. 

Corresponding  to  Equation  2.2,  we  can  define  a  near  zone  field  as 


£V(r,M)  =  n^£(r,0,*).  (2.4) 


Note  that  Ex  is  a  function  of  r  as  well  as  6  and  <£,  since  for  finite  r,  the 
radial  dependence  of  E  ia  not  exactly  e"*' /r.  However,  as  r  gets  large,  Ex 
will  become  independent  of  r  and  approach  By.  We  now  define  the  near 
zone  gain  as 

C»jv(r,0.tf)  -  iQp  *  (2-5) 


The  advantage  of  the  above  definition  of  near  zone  gain  is  that  it  allows 
Gft  to  be  directly  compared  to  G>  without  the  need  for  a  radial  dependent 
scale  factor.  Also,  as  r  gets  large  G x  will  approach  Gy. 

If  in  Equation  2.3  we  insert  the  9  component  of  the  electric  field,  then 
we  obtain  the  9  polarized  gain,  If  we  insert  the  $  polarization,  then  we  get 
(lie  d  polarized  gain.  Similarly,  the  near  zone  gain  of  Equation  2.5  has  r, 
9,  and  ^  components. 

Next  consider  the  definition  of  RCS  in  the  near  and  far  zone.  We  will 
still  denoie  the  radiated  field  of  the  target,  b«l  we  will  now 


assume  that  the  target  currents  were  induced  by  a  plane  wa  •  of  magnitude 
E,  incident  from  the  direction  (#,,<?,).  The  far  zone  bistatic  echo  area  or 
RCS  is  given  by 


=  4* 


|£,!: 


(2.6) 
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If  (0,  <f>)  =  (${,<$>,),  then  <rF  is  referred  to  as  the  backscatter  RCS.  If  (0,  $)  = 
(ir  -  0,,-jt  +  <pt),  then  crF  is  referred  to  as  the  forward  RCS.  Similarly  we  will 
define  the  near  zone  bistatic  RCS  as 


<7w(r,0,d>)  =  4rr 


\Et\> 


(2.7) 


< TF  has  four  components,  corresponding  to  the  possible  polarizations  of 
the  incident  and  far  zone  scattered  plane  waves.  These  are: 

6  incident  and  6  scattered 

0  incident  and  0  scattered 

0  incident  and  0  scattered 

0  incident  and  0  scattered. 

<r(V  has  the  above  four  components  plus 

B  incident  and  f  scattered 

0  incident  and  r  scattered 

since  the  near  zone  scattered  field  can  have  a  radial  component.  ESP  IV 
will  output  the  magnitude  and  phase  of  the  near  or  far  zone  gain  and  RCS. 
By  the  phase  of  the  near  or  far  zone  gaiu  or  RCS,  we  mean  the  phase  of 
Eft  or  E?, respectively. 

The  computation  of  near  zone  fields  will  be  illustrated  with  the  four 
element  dipole  array  shown  in  Figure  2.1.  Here  we  have  four  A/2  dipoles 
located  A/4  in  front  of  a  perfectly  conducting  1.2A  x  Q.75A  reflecting  plate. 
A  three  view  plot  of  the  wire  and  plate  geometry,  supplied  by  tht  ESP  code, 
is  shown  in  Figure  2.2.  Here  it  can  lx*  seen  that  for  the  MM  solution  each 
dipole  is  split  into  4  equal  segments,  w'hich  results  in  3  piecewise  sinusoidal 
modes  on  each  dipole.  The  plate  is  modelled  by  38  surface  patch  modes. 
The  HipnW  ere  perfectly  eondnetinjj  with  a  radius  of  O.flfll  V.  F.ft»*h  di)»»*|e  ic 
fed  by  an  in  phase  one  volt  generator,  which  in  the  far  zone  will  produce  a 
beam  maximum  on  the  -far  axis,  that  is  at  0  =  90°  and  0  =  0.  We  will  now 
show  that  as  the  radius  of  the  near  zone  gain  patterns  for  the  dipole  array 
increases,  the  near  tone  gain  patterns  approach  the  far  zone  gain  patterns. 
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I.EX  X  0,75 X  PLATE 


/ 


Figure  2.1:  Top  View  of  four  A/2  dipoles  located  A/4  iu  front  of  a 
1.2A  x  0.75A  reflecting  plate. 

F-ir  the  dipole  array,  Figures  2.3, 2.4  show  the  far  cone  gain  iu  the  az¬ 
imuth  plane  fi  =  90°  and  in  the  elevation  plane  $  »  0.  Only  the  B  polarita- 
tion  is  shown,  since  in  these  planes  there  is  no  polarisation.  By  eoutrasi 
Figures  2.5-2.9  show  the  $  polarised  near  rone  gain  in  the  asinmth  plane 
and  for  r  =  0.7,  1,  2,  5,  and  10  A.  Figures  2.10*2.14  and  2.15*2.19  show  the 
&  and  r,  polarised  new  tone  gain  in  the  elevation  plane  and  for  the  same  r 
values.  Note  that  all  of  these  pattern  plots  are  supplied  by  the  ESP  code. 

Comparing  Figures  2.3  and  2.5  shows  that  the  near  roue  azimuth  plane 
pattern  at  r  «  0.7A  is  considerably  different  from  the  far  tone  pattern. 
In  particular,  the  peak  gain  at  r  =  Q.7A  is  7.2  dB,  as  compared  to  9.3 
dB  in  the  far  rone.  Also,  the  beamwidth  at  r  =  0.7A  is  much  wider  than 
the  beamwidth  in  the  far  tone.  There  arc  two  main  reasons  for  these 
differences.  First,  the  uniform  phase  taper  will  produce  radiation  from 
each  dipole  which  adds  iu  phase  at  $  =  90°  and  ^  »  0  in  the  far  tone  tone, 
but  not  M  r  a  OTA.  Thus,  the  peak  near  tone  gain  is  less  than  the  peak 
far  tone  gain.  Second,  in  the  near  tone  the  end  elements  in  the  dipole  array 
are  much  closer  to  the  field  points  near  $  =  db90°  than  the  center  elements. 
This  increases  the  field  intensity  near  d  =  &90°  and  thus  increases  the  near 
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GEOMETRY  FOR  4  ELEMENT  DIPOLE  ARRAY 


12  WIRE  MODES 
38  PLRTE  MODES 
0  RTTflCH.  MODES 
50  TOTAL  MODES 
[SCRLE1  =  0.440  A 


Z  AXIS  VIEW 


X  AXIS  VIEW  Y  AXIS  VIEW 


Figure  2.2:  A  three  view  sketch  of  the  four  element  dipole  array  and  re¬ 
flecting  plate  of  Figure  2.1. 
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zone  beamwidth.  Comparing  Figure  2.3  to  Figures  2.S-2.9  shows  that  as  r 
increases,  the  near  zone  azimuth  patterns  approaches  the  far  zone  pattern. 

Figures  2.10-2.14  show  the  0  polarized  near  zone  elevation  plane  patterns 
at  r  =  0.7, 1, 2,  5,  and  10  A.  Comparing  Figures  2.4  and  2.10  shows  that  the 
near  zone  elevation  plane  pattern  at  r  =  0.7A  is  considerably  different  than 
the  far  zone  pattern.  However,  comparing  Figure  2.4  to  2.10-2.14  shows 
chat  as  r  increases  the  near  zone  elevation  plane  patterns  approach  the  far 
zone  pattern. 

In  the  far  zone  the  f  component  of  aJl  fields  is  zero,  but  in  the  near  zone 
it  is  not.  For  the  dipole  array  Figures  2.15-2.19  show  the  f  polarized  gain 
in  the  near  zone  and  for  r  —  0.7,  I,  2,  5,  and  10  A.  Note  that  as  r  increases, 
t  he  level  of  the  ?  polarized  gain  decreases  from  -6.4  dB  at  r  =  0.7A  to  -22.6 
dB  at  r  —  10A. 

Next  we  will  consider  the  RCS  of  th**  dipole  array  shown  in  Figures 
2.1  and  2.2.  For  a  $  polarized  incident  and  scattered  plane  wave,  Figure 
2.20  shows  the  far  zone  backscatter  RCS  pattern  in  the  azimuth  plane 
0  -  90°.  Figures  2.21-2.25  show  the  corresponding  near  zone  RCS  patterns 
for  r  =  0.7,  1,  2,  5,  and  10  A.  Although  the  RCS  is  shown  in  units  of 
dB  over  a  square  meter,  in  this  case  it  is  equivalent  to  dB  over  a  square 
wavelength,  since  the  computations  were  made  at  a  frequency  of  300  Mhz 
(A  =  1  meter).  Comparing  Figures  2.20  and  2.21  shows  that  the  RCS  at 
r  =  G.7A  is  considerably  different  from  the  far  zone  RCS.  However,  as  r 
increases,  the  near  zone  RCS  does  approach  the  far  zone  RCS. 

In  summary,  the  capability  to  compute  near  zone  fields  has  been  added 
to  the  ESP  code.  Other  than  numerical  integrations,  no  approximations 
are  involved  in  the  computation  of  the  near  zone  fields,  which  are  not  also 
made  in  the  computation  of  the  far  zone  fields.  The  above  data  suggests 
that  we  have  properly  modified  ESP  to  compute  near  zone  fields.  The 
main  limitations  are  that  if  the  field  point  gets  too  close  to  the  wire/plate 
geometry,  artificial  slope  discontinuities  in  the  MM  expansion  modes  will 
cause  anomalous  behavior  in  the  computed  near  zone  fields.  For  this  reason 
the  field  point  should  not  be  closer  than  about  a  tenth  of  a  wavelength  from 
the  wirc/plate  geometry. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 


GRIN 

POLAR  I ZRT I  ON:  0 

flZIM.  PLANE:  0  =90°  R  =  00 

MAXIMUM  =  9.300  DB  10  DB/DIV. 


Figure  2.3:  0  polarized  far  zone  gain  for  the  dipole  array  in  the  azimuth 
plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLARIZATION:  6 

ELEV.  PLANE : 4  =  0°  R  =  « 

MAXIMUM  =  9.300  DB  10  DB/DIV. 


Figure  2  4:  0  polarized  far  zone  gain  for  the  dipole  array  in  the  elevation 
plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLARIZATION:  0 

AZIM.  PLANE: 0  =  90°  R=  0.700  M 
MAXIMUM  =  7.200  OB  10  DB/DIV. 

-h  $ 


Figure  2.5:  9  polarized  near  zone  gain  at  r  =  0.7A  for  the  dipole  array  in 
the  azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GAIN 

POLARIZATION:  6 

AZIM.  PLANE: 0  =  90“  R=  1.000  M 
MAXIMUM  =  8.300  DB  10  DB/DIV. 

-U.  (j) 


Figure  2.6:  B  polarized  near  zone  gain  at  r  =  1A  for  the  dipole  array  in  the 
azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLARIZATION:  0 

flZIM.  PLANE: ®  =  90”  R=  2.000  M 
MAXIMUM  =  9.100  DB  10  DB/DIV. 

+-  $ 


Figure  2.7:  6  polarized  near  zone  gain  at  r  =  2A  for  the  dipole  array  in  the 
azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GAIN 

POLARIZATION:  6 

flZIM.  PLANE: 0  =  90°  R=  5.000  N 
MAXIMUM  =  9.300  DB  10  DB/DIV. 


Figure  2.8:  6  polarized  near  zone  gain  at  r  =  5A  for  the  dipole  array  in  the 
azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  hHZ. 

GRIN 

POLAR  I ZRT I  ON:  6 

AZIN.  PLANE: @  =  90°  R=  10.000  M 
MAXIMUh  =  9.300  DB  10  DB/DIV. 

-f-  $ 


* 

Figure  2.9:  0  polarized  near  zone  gain  at  r  =  10A  for  the  dipole  array  in 
the  azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GAIN 

POLARIZATION:  0 

ELEV .  PLANE: $  =  0°  R=  0.700  M 

MAXIMUM  =  7.200  DB  10  OB/DIV. 


—  0 


Figure  2.10:  6  polarized  near  zone  gain  at  r  —  0.7A  for  the  dipole  array  in 
the  elevatiou  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLAR  I ZRT I  ON:  % 

ELEV .  PLANE :  4  =  0°  R=  1.000  M 

MAXIMUM  =  8.300  D8  10  OB/DI V. 


—  0 


Figure  2-11:  0  polarised  near  Zone  gain  at  r  =  1A  for  the  dipole  array  in 
the  elevation  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLARIZATION:  0 

ELEV .  PLANE : $  =  0°  R=  2.000  M 

MAXIMUM  =  9.100  OB  10  DB/'DIV. 


Figure  2.12:  0  polarized  near  zone  gain  at  r  =  2A  for  the  dipole  array  in 
the  elevation  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

PQLfiRIZRTION:  6 

ELEV.  PLANE : 4  =  0°  R=  5.000  M 

MAXIMUM  =  9.300  DB  10  DB/DIV. 


Figure  2.13:  &  polarized  aear  tone  gain  at  r  ^  5A  for  the  dipole  array  ic 
the  elevation  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLAR  I ZRT I  ON:  6 

ELEV.  PLANE: *  =  0°  R=  10.000  M 

MAXIMUM  =  9.300  DB  10  OB/DIV. 


Figure  2.14:  $  polarised  near  sone  gain  at  r  »  10A  for  the  dipole  array  in 
the  elevation  plane- 
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4  ELEMENT  DIPOLE  ARRAY  ' 

FREQUENCY  =  300.00  MHZ. 

GAIN  .  .  ■ 

POLARIZATION:  '  R  .  / 

ELEV .  PLANE:  4>  =  0°  "R=  0.700  M 

MAX I MUM  =  -6.400  DB  10  DB/DIV. 


Figure  2.15:  f  polarized  near  zone  gain  at  r  —  0.7A  for  the  dipole  array  in 
the  elevation  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GAIN 

POLARIZATION:  R 

ELEV .  PLANE : $  =  0°  R=  1.000  M 

MAXIMUM  =  -6.300  DB  10  DB/DIV. 


—  0 


Figure  2.16:  r  polarized  near  zone  gain  at  r  s=  1A  for  the  dipole  array  in 
the  elevation  plane. 


4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLARIZATION:  R 

ELEV.  PLANE: $  =  0°  R=  2.000  M 

MAXIMUM  =  -9.600  DB  10  DB/DIV. 


—  0 


Figure  2.17:  r  polarized  near  zone  gain  at  r  2A  for  the  dipole  array  in 
the  elevation  plant*. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300 . 00  MHZ'. 

GAIN 

POLARIZATION:  -  R 

ELEV. 'PLANE:  4  =  0°  R=  5.000  ‘N 

MAXIMUM  =-16.700  DB  10  DB/DIV. 


0 


Figure  2.18:  r  polarized  near  zone  gain  at  r  =  5X  for  the  dipole  array  in 
the  elevation  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 

GRIN 

POLARIZATION:  R 

ELEV .  PLANE :  4  =  0°  R=  10.000  M 

MAXIMUM  =-22.600  DB  10  DB/DIV. 

4"*‘  0 


Figure  2.19:  f  polarized  near  zone  gain  at  r  =  10A  for  the  dipole  array  in 
the  elevation  plane. 
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4  ELEMENT  DIPOLE  ARRAY 

FREQUENCY  =  300.00  MHZ. 
BRCKSCHTTER 

POLARIZATION:  6-IN  0-OUT 
flZIH.  PLANE:®  -  90°  R  =  » 
MAXIMUM  =  11.500  DB/M2  10  DB/DIV. 


Figure  2.20:  The  §  polarized  far  zone  backscatter  RCS  for  the  dipole  array 
in  the  azimuth  plane. 
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4  ELEMENT  BIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 
BFICKSCfiTTER 

POLAR I ZRT I ON:  6-IN  0-OUT 
AZIM.  PLANE:®  =  90°  R-  0.700  M 
MAXIMUM  =  9.500  DB/M2  10  DB/DIV. 

-u  $ 


Figure  2.21:  The  $  polarized  near  zone  backscatter  RCS  for  the  dipole  array 
at  r  =  0.7 A  in  the  azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 
BflCKSCRTTER 

POLARIZATION:  0-IN  0-OUT 
R7IN  PLANE:®  =  90°  R=  1.000  M 
MAXIMUM  =  10.600  DB/M2  10  DB/OIV. 

+*  $ 


Figure  2.22:  The  §  polarized  near  zone  backscatter  RCS  for  the  dipole  array 
at  v  =  A  iu  the  azimuth  plane. 


4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 
8RCKSCRTTER 

POLARIZATION:  0-IN  0-OUT 
flZIM.  PLANE : 0  =  90°  R=  2.000  M 
MAXIMUM  =  11.500  DB/M2  10  DB/DIV. 


Figure  2.23:  The  0  polarized  near  zone  backscatter  RCS  for  the  dipole  array 
at  r  =  2A  in  the  azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  MHZ. 
BflCKSCATTER 

POLARIZATION:  0-IN  0-OUT 
nzili.  PLANE : 0  =  90°  R=  5.000  M 
MAXIMUM  =  11.600  DB/M2  10  DB/DIV. 

+-  $ 


Figure  2.24:  The  9  polarized  near  zone  backscatter  RCS  for  the  dipole  array 
«  a=  5A  in  the  azimuth  plane. 
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4  ELEMENT  DIPOLE  ARRAY 
FREQUENCY  =  300.00  NHZ. 
BflCKSCRTTER 

POLRRIZflTION:  0-IN  0-OUT 
flZIM.  PLANE: 0  =  90°  R=  10.000  M 

nnxinuii  =  11.600  db/m2  10  db/div. 


Figure  2.25:  The  $  polarised  near  zone  bacfcscatter  RCS  for  the  dipole  array 
at  r  —  10A  in  the  azimuth  plane. 
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Figure  2.26:  Side  view  of  &  thin  dielectric  plate. 

2,2  Thin  Dielectric  Plates 

Previous  versions  of  the  ESP  have  been  limited  to  treating  plates  of  zero 
thickness  and  perfect  conductivity.  In  this  section  we  will  describe  the  use 
of  the  sheet  impedance  approximation  which  will  permit  version  IV  of  BSP 
to  treat  electrically  thin  dielectric  plates. 

Figure  2.26  shows  the  side  view  of  a  dielectric  plate  of  thickness  T  and 
relative  dielectric  constant  tr  =  </*o  immersed  in  free  space  with  permit¬ 
tivity  €0  %  fi.85  x  IQ"13  F/m.  If  the  slab  is  lossy,  then  er  will  be  complex. 
The  sheet  impedance  approximation  is  valid  provided: 

1,  The  dielectric  plate  is  sufficiently  thin  that  the  electric  field  is  essen¬ 
tially  uniform  throughout  its  thickness.  This  will  be  the  case  provided 
Iky^TI  <  1,  where  k  is  the  free  space  wavenumber. 

2.  The  dominant  polarization  of  the  electric  field  in  the  dielectric  slab  is 
parallel  to  the  broad  surfaces  of  the  slab.  This  assumption  is  usually 
best  satisfied  when  the  incident  field  is  polarized  parallel  to  the  broad 
surfaces  of  the  slab  or  if  |c,|  >  1. 
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If  the  above  assumptions  are  satisfied,  then  the  thin  dielectric  plate  can 
be  modeled  as  a  zero  thickness  plate  with  sheet  impedance  [4]-[8] 

Z‘  =  ,*«,(<.  -  1)T 

The  sheet  impedance  of  a  multilayer  slab  is  simply  the  parallel  combination 
of  the  sheet  impedance  of  each  layer  (8).  For  example,  if  we  have  a  two  layer 
slab,  and  if  Zt\  and  Z,i  denote  the  sheet  impedance  of  the  individual  layers, 
then  the  sheet  impedance  of  the  two  layer  slab  is 


Z$\  4-  Zt j 


(2.9) 


From  Equation  2.8  it  can  be  seen  that  if  the  dielectric  slab  is  lossless,  then 
Zt  is  pure  imaginary.  If  the  slab  is  very  lossy,  such  that  lm(er)  >  Re(«,), 
then  Equation  2.8  becomes 


(2,10) 


where  o  is  the  conductivity  of  the  material.  A  material  with  pure  real  sheet 
impedance,  as  given  by  Equation  2.10,  is  referred  to  as  a  resistive  sheet. 

Probably  the  main  advantage  of  the  sheet  impedance  approximation  is 
that  it  is  relatively  simple  (as  compared  to  more  exact  models)  to  imple¬ 
ment.  In  particular,  assume  one  has  a  surface  patch  MM  computer  code 
for  perfectly  conducting  plates,  such  as  the  ESP  code.  Then  to  modify  this 
code  to  treat  plates  with  a  sheet  impedance,  it  is  only  necessary  to  add  the 


to  the  mn  term  in  the  MM  impedance  matrix  for  the  perfectly  conducting 
plates.  In  Equation  2.11,  JR,  is  the  MM  basis  function,  and  the  in¬ 
tegration  is  over  that  portion  of  the  surface  where  Jw,  Ju>  and  Zt  are  all 
non-zero.  Not#-  that  the  integration  in  Equation  2.11  is  reasonably  simple, 
even  if  Z,  varies  over  the  surface,  if  Z,  —  0,  then  =  0.  A  sheet 

impedance  plate  is  thus  a  generalization  of  the  perfectly  conducting  plate, 
with  the  perfectly  conducting  plate  being  the  special  case  Zt  —  0.  ESP 
Version  IV  has  beer,  set  up  to  allow  the  user  to  specify  a  different  Zt  for 
each  plate.  Thus,  ESP  IV  will  be  able  to  treat  combinations  of  perfectly 
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conducting  and  thin  dielectric  plates.  Also,  ESP  IV  will  be  able  to  compute 
the  near  or  far  zone  fields  of  a  surface  impedance  plate. 

In  our  present  implementation  of  the  MM  solution  for  a  sheet  impedance 
plate,  we  use  the  same  current  expansion  or  basis  functions  as  for  perfectly 
conducting  plates.  On  a  perfectly  conducting  plate  the  normal  component 
of  the  surface  current  density  goes  to  zero  at  the  edges  of  the  plate,  and  our 
basis  functions  enforce  this  condition.  However,  on  a  sheet  impedance  plate, 
the  normal  component  of  the  surface  current  density  need  not  vanish.  Thus, 
the  basis  functions  presently  being  used  for  sheet  impedance  plates  are  not 
optimum.  The  solution  to  this  problem  is  to  augment  our  present  expansion 
modes  with  a  series  of  edge  modes  which  allow  for  a  finite  normal  component 
of  current  at  the  edges  of  the  sheet  impedance  plates.  As  described  below, 
this  has  been  done  for  a  single  plate  (see  Figures  2.36,2.37  below),  but  a 
general  technique  not  not  been  implemented  in  the  ESP  code. 

The  first  example  of  the  use  of  the  sheet  impedance  model  for  a  thin 
dielectric  plate  is  shown  in  the  insert  in  Figure  2.27.  Here  we  have  a  0.813 
A  square  dielectric  plate  of  thickness  T  =  0.051  A  and  relative  dielectric 
constant  tr  ~  4,  Using  Equation  2.8,  this  dielectric  plate  has  a  sheet 
impedance  of  Z,  =  -;302  fl.  Figure  2.27  shows  a  comparison  of  the  mea¬ 
sured  (9j  and  computed  RCS  in  the  elevation  plane  and  for  4>  polarization. 
The  same  data  is  shown  for  the  9  polarization  in  Figure  2.28.  Although 
there  is  overall  agreement  between  computed  and  measured  results  in  these 
figures,  there  are  also  noticeable  differences.  One  reason  for  the  differences 
is  that  for  this  dielectric  plate 


kyfCT  =  0.64 


which  is  somewhat  large  to  expect  highly  accurate  results  from  the  sheet 
impedance  model.  Also,  for  the  B  polarized  RCS  of  Figure  2.28,  near  grating 
incidence  (8,  =s  0,  180°)  the  incident  field  is  polarized  perpendicular  to 

the  broad  surfaces  of  the  plate.  Thus,  in  this  case  the  dominant  electric 
field  in  the  dielectric  will  be  perj>cndicular  to  the  broad  surfaces  of  the 
slab.  As  described  above,  this  violates  the  second  condition  for  the  sheet 
impedance  approximation.  We  note  that  the  presence  of  the  elertric  field 
perpendicular  to  the  broad  surface  of  the  dielectric  plate  con  be  accounted 
for  by  a  magnetic  sheet  current,  however,  the  complexity  of  the  resultiug 
solution  is  greatly  increased  (10). 
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Figure  2.27:  The  RCS  of  a  0.813  A  square  dielectric  plate  at  300  Mhz.  and 
for  <p  polarization 
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The  second  example  of  the  use  of  the  sheet  impedance  model  is  illus¬ 
trated  in  the  insert  in  Figure  2.29.  Here  we  have  a  dielectric  cylindrical  shell 
of  height  3.73  cm,  outer  radius  1.27  cm,  shell  thickness  T  =  0.254  cm,  and 
relative  dielectric  constant  er  —  2.54.  Using  Equation  2.8,  at  /  =  6.03  Ghz 
this  dielectric  shell  is  equivalent  to  the  sheet  impedance  Z,  —  —]763  Cl. 
For  the  purpose  of  the  MM  model,  the  circular  cylindrical  dielectric  shell 
is  modelled  by  a  zero  thickness  octagonal  cylinder  of  radius  1.14  cm  and 
surface  impedance  — ;763  Cl.  56  surface  patch  modes  are  used  to  model 
the  current  on  the  cylinder. 

Figures  2.29  and  2.30  show  the  dielectric  cylinder  elevation  plane  RCS 
patterns  for  the  0  and  <f>  polarizations.  In  these  curves  the  solid  line  is 
measurements  by  Van  Doeren  [11],  the  dots  are  computations  by  Van  Do- 
eren,  and  the  x’*  are  our  sheet  impedance  model  computations.  For  this 
dielectric  shell 

ky/TtT  =  0.32. 

The  agreement  between  our  and  Van  Doeren ’s  results  is  only  fair. 

As  a  final  example,  Figure  2.31  shows  a  dielectric  strip  of  length  3.048 
m,  width  0.1829  m,  and  thickness  T  =  0.0001829  m,  and  with  er  =  4.  The 
dielectric  strip  is  excited  by  a  unit  amplitude  zori  polarized  plane  wave 
incident  from  the  +y  axis,  i.e., 

E[  =  El  =  ejkv  V/m. 

At  /  ~  1  GHz,  this  thin  dielectric  strip  has  the  equivalent  sheet  impedance 
Z,  =  ~;32800  Cl. .  For  a  z  polarized  incident  and  scattered  wave  Figure 
2.32  and  2.33  show  the  magnitude  and  phase  of  the  near  zone  scattered 
electric  field  along  the  y  axis  (Note:  If?]  in  dB  =  201og|.5?|  in  V/m),  The 
solid  line  is  a  computation  from  Gupta  [12]  and  the  dashed  line  is  the  sheet 
.impedance  model.  Figures  2.34  and  2.35  show  the  same  data  for  the  x 
polarization.  In  brief  Gupta  uses  a  volume  polarization  current  model  for 
the  dielectric  strip,  and  in  this  way  is  more  exact  than  our  sheet  impedance 
model.  However,  he  also  uses  asymptotic  and  far  zone  approximations 
in  finding  the  fields  of  the  currents,  and  in  this  wav  is  less  exact  than  onr 
model.  These  Figures  show  almost  perfect  agreement  for  the  z  polarization, 
while  the  magnitude  of  the  x  polarized  scattered  field  is  about  1  to  2  dB 
below  Gupta’s  computation. 

The  above  three  examples  illustrate  that  our  implementation  of  the 
sheet  impedance  model  is  capable  of  producing  results  within  a  few  dB 
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Figure  2.29:  The  backscatter  RCS  of  r  dielectric  cylindrical  shell  for  0 
polarization  of  the  incident  electric  held. 
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Figure  2.30:  The  backscatter  RCS  of  a  dielectric  cyiiudrical  shell  for  ^ 
polarization  of  the  incident  electric  field. 
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Figure  2.31:  Side  view  of  an  electrically  thin  dielectric  plate. 


of  the  correct  result.  Although  these  results  may  be  accurate  enough  for 
engineering  accuracy,  some  improvement  is  possible.  In  particular  consider 
Figure  2.34,  where  our  results  are  a  couple  of  dB  below  Gupta’s  computa¬ 
tions.  For  this  dielectric  strip 

ky/erT  =  0.0011, 

and  the  incident  field  is  polarized  parallel  to  the  broad  face  of  the  strip. 
Thus  we  would  expect  the  sheet  impedance  approximation  to  be  very  ac¬ 
curate,  and  would  not  expect  the  2  dB  difference  evident  in  Figure  2.34. 
Thus,  as  described  above,  we  expect  that  a  substantial  portion  of  the  error 
in  Figure  2.34  is  due  to  the  fact  that  our  expansion  for  the  surface  current 
density  on  the  sheet  impedance  plates  vanishes  at  the  edges,  while  the  ac¬ 
tual  current  does  not.  To  test  this  hypothesis,  we  added  a  series  of  edge 
modes  which  allowed  for  finite  surface  current  density  at  the  edge  of  the 
sheet  impedance  plate.  Figures  2.36,2.37  repeat  the  same  data  as  in  Figures 
2.34,2.35,  except  the  computation  included  these  edge  modes.  Comparing 
Figures  2.34  and  2.36  shows  that  adding  the  edge  modes  considerably  im¬ 
proved  the  agreement  between  our  and  Gupta’s  results. 

In  summary  we  have  implemented  in  the  BSP  code  a  simple  sheet  for 
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thin  dielectric  plates.  The  simple  sheet  impedance  model  is  capable  of 
predicting  near  or  far  zone  fields  to  within  a  few  dB  of  the  actual  values. 
The  accuracy  of  our  model  could  be  improved  through  the  addition  of  the 
above  described  edge  modes. 
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Figure  2,32:  The  magnitude  (V/m)  of  the  near  zone  Ez  field,  on  the  line 
( x  =  0 ,y,z  =  0),  for  the  dielectric  strip  of  Figure  2.31  caused  by  a  normally 
incident  i  polarized  plane  wave. 
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Figure  2.33:  The  phase  of  the  near  zone  Et  field,  on  the  line 

(*  =  V> z  ~  °)>  the  dielectric  strip  of  Figure  2.31  caused  by  a  normally 
incident  z  polarized  plane  wave. 
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Figure  2-34:  The  magnitude  (V/m)  of  the  near  sons  Et  field,  on  the  line 
(z  =s  0,  y,  z  =a  0),  for  the  dielectric  strip  of  Figure  2.31  caused  by  a  normally 
incident  z  polarized  plane  wave. 
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Figure  2.35:  The  phase  of  the  near  zone  Ex  field,  on  the  line 
(x  =  0,  y, x  =  0),  for  the  dielectric  strip  of  Figure  2.31  caused  by  a  normally 
incident  i  polarised  plane  wave. 
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Figure  2.36:  The  magnitude  (V/m)  of  the  near  zone  Es  field,  on  the  line 
(*  =  0,y,  2  =  0),  for  the  dielectric  strip  of  Figure  2.31  roused  by  a  normally 
incident  i  polarized  plane  wave.  Data  computed  with  edge  modes. 
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Chapter  3 
Summary 


This  report  lias  described  two  recent  modifications  to  the  ESP  code.  ESP 
is  a  user  oriented  code  for  the  analysis  of  the  electromagnetic  radiation  and 
scattering  from  geometries  consisting  of  interconnections  of  ihir  wires  and 
plates.  The  two  modifications  or  improvements  are: 

1.  the  ability  to  compute  near  zone  radiation  and  scatter  ng  patterns 
(previously  ESP  could  only  compute  far  zone  patterns), 

2.  the  ability  to  analyze  thin  dielectric  plates  through  a  shee  impedance 
model  (previously  ESP  could  only  analyze  perfectly  conducting  plates). 

These  two  features  should  be  pari  of  ihe  next  version  (IV)  of  the  ESP  code. 
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MISSION 

°f 

Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (CSI)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  C*/  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


